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Abstract: We investigate an oblivious routing scheme, amenable to distributed computation and 
resilient to graph changes, based on electrical flow. Our main technical contribution is a new round- 
ing method which we use to obtain a bound on the C\ — > l\ operator norm of the inverse graph 
Laplacian. We show how this norm reflects both latency and congestion of electric routing. 



1 Introduction 

Overview We address a vision of the Internet where every participant exchanges messages with their 
direct friends and no one else. Yet such an Internet should be able to support reliable and efficient 
routing to remote locations identified by unchanging names in the presence of an ever changing graph 
of connectivity. 

Modestly to this end, this paper investigates the properties of routing along the electric flow in a 
graph {electric routing for short) for intended use in such distributed systems, whose topology changes 
over time. We focus on the class of expanding graphs which, we believe, gives a good trade-off between 
applicability and precision in modeling a large class of real-world networks. We address distributed 
representation and computation. As a measure of performanace, we show that electric routing, being 
an oblivious routing scheme, achieves minimal maximum edge congestion (as compared to a demand- 
dependent optimal scheme). Furthermore, we show that electric routing continues to work (on average) 
in the presence of large edge failures ocurring after the routing scheme has been computed, which attests 
to its applicability in changing environments. We now proceed to a formal definition of oblivious routing 
and statement of our results. 
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Oblivious routing The object of interest is a graph G = (V, E) (with V = [n] and \E\ = m) undirected, 
positively edge-weighted by w Ufi > 0, and not necessarily simple. The intention is that higher w UtV 
signifies stronger connectivity between u and v; in particular, w UjV = indicates the absence of edge 
(u, v). For analysis purposes, we fix an arbitrary orientation "— >" on the edges (u, v) of G, i.e. if (u, v) 
is an edge then exactly one of u — > v or v — > u holds. Two important operators are associated to every 
G. The discrete gradient operator B e R ExV , sending functions on V to functions on the undirected 
edge set E, is defined 2&x\ UJ0 ) B := Xu ~Xv if u -> v, mdx* {uv) B := %v ~ Xu otherwise, where x y is the 
Kronecker delta function with mass on y. For e e E, we use the shorthand B e :- ix e B)* . The discrete 
divergence operator is defined as B*. 

A (single-commodity) demand of amount a > between s e V and t e V is defined as the vector 
d - a(xs - Xt) <= R V - A (single-commodity) flow on G is defined as a vector / e R £ , so that /(„ „) 
equals the flow value from u towards v if u — > u, and the negative of this value otherwise. We also use 
the notation := f( Ufi) if m — > v, and /„^„ := -f( u , v ) otherwise. We say that flow / routes demand d 
if B*f = d. This is a linear algebraic way of encoding the fact that / is an (s, ?)-flow of amount a. A 
multi-commodity demand, also called a demand set, is a matrix whose columns constitute the individual 
demands' vectors. It is given as the direct product ® T d T of its columns. A multi-commodity flow is 
represented as a matrix © T / T , given as a direct product of its columns, the single-commodity flows. For 
clarity, we write f T , e for (f T ) e . The flow © T / T routes the demand set ® T d T if B*f T = d T , for all t, or in 
matrix notation B*(® T f T ) = ® T d T . The congestion \\ • \\c of a multi-commodity flow measures the load of 
the most-loaded edge, relative to its capacity. It is given by 

||© T / T || G :=maxV \f T , e /w e \ = \\(® T f T yW-%^, where \\A\\^ : =sup^. (1.1) 
e V x *o IWIi 



An oblivious routing scheme is a (not necessarily linear) function R : K v — > R £ which has the 
property that R(d) routes d when d is a valid single-commodity demand (according to our definition). 
We extend R to a function over demand sets by defining R(® T d T ) := <S T R(d T ). This says that each demand 
in a set is routed independently of the others by its corresponding R-How. We measure the "goodness" 
of an oblivious routing scheme by the maximum traffic that it incurs on an edge (relative to its capacity) 
compared to that of the optimal (demand-dependent) routing. This is captured by the competitive ratio 
t]r of the routing scheme R, defined 

\\R{®rd T )\\ G 

% :=sup sup — — . (1.2) 

® T d T ® T f T II ©T Jt\\G 

B'(® T f T )=® T d T 

Let £ denote the (yet undefined) function corresponding to electric routing. Our main theorem states: 
THEOREM 1.1 — For every undirected graph G with unit capacity edges, maximum degree d max 

C | 

and vertex expansion a :- minjcv ~[j|j^fjj' one ^ as ^ ^ {^^ n f) ' { a ^ n idt^-a ) • •^ u ' s JS tight up to 
a factor of 0(ln In n). 
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The competitive ratio in Theorem 1 . 1 is best achievable for any oblivious routing scheme up to a 
factor of 0(lnln?i) due to a lower bound for expanders, i.e. the case a = 0(1), given in [7]. Theorem 1.1 
can be extended to other definitions of graph expansion, weighted and unbounded-degree graphs. We 
omit these extensions for brevity. We also give an unconditional, albeit much worse, bound on rjf. 

THEOREM 1.2 — For every unweighted graph on m edges, electrical routing has rj£ < 0(m l l 2 ). 
Furthermore, there are families of graphs with corresponding demand sets for which 77 £ - Q.{m 1 ^ 2 ). 

Electric routing Let W = diag(. . . , w e , . . . ) e R ExE be the edge weights matrix. We appeal to a known 
connection between graph Laplacians and electric current [5, 13]. Graph edges are viewed as wires of 
resistance w~ l and vertices are viewed as connection points. If <p e R y is a vector of vertex potentials 
then, by Ohm's law, the electric flow (over the edge set) is given by / = WB<p and the corresponding 
demand is B*f = Lip where the (un-normalized) Laplacian L is defined as L = B*WB. Central to the 
present work will be the vertex potentials that induce a desired (s, ?)-flow, given by - V(x s - Xt), 
where iJ is the pseudo-inverse of L. Thus, the electric flow corresponding to the demand pair (s, t) is 
WB(p^ = WBt}(Xs - Xt)- We define the electric routing operator as 



The vector E(x„ ~Xt) 6 R £ encodes a unit flow from s to t supported on G, where the flow along an edge 
(u, v) is given by \st, uvj := £(^ 4 - Xt)u->v = (fl?'^ ~ $'^) w u,v X (O ur convention is that current flows 
towards lower potential.) When routing an indivisible message (an IP packet e.g.), we can view the unit 
flow E(Xs ~Xt) as a distribution over (s, £)-paths defined recursively as follows: Start at s. At any vertex 
u, forward the message along an edge with positive flow, with probability proportional to the edge flow 
value. Stop when t is reached. This rule defines the electric walk from s to t. It is immediate that the 
flow value over an edge (u, v) equals the probability that the electric walk traverses that edge. 

Let "~" denote the vertex adjacency relation of G. In order to make a (divisible or indivisible) 
forwarding decision, a vertex u must be able to compute \st, uv\ for all neighbors v ~ u and all pairs 



Representation In order to compute \st, uv\ (for all s, t e V and all v ~ u) at u, it suffices that u stores 
the vector tp^ := L*x w , for all w e {w : w ~ u] U {u}. This is apparent from writing 



where we have (crucially) used the fact that is symmetric. The vectors tp^ stored at u comprise the 
(routing) table of u, which consists of deg(w) • n real numbers. Thus the per- vertex table sizes of our 
scheme grow linearly with the vertex degree - a property we call fair representation. It seems that fair 
representation is key for routing in heterogenous sytems consisting of devices with varying capabilities. 

'The bilinear form fst, uvj = Xs,iBL*B"Xu,t acts like a "representation" of G, hence the custom bracket notation. 



E(d) = WBLU 



(1.3) 




lst,uv\ = iXu -Xv)L f (x s -Xt) = (<p [u] - <p lv] )*(x s -Xt), 



(1.4) 
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Equation (1.4), written as [jf, hu] - (Xs~Xt)*(<P <P )> shows that in order to compute [sf, hu] atw, 
it suffices to know the indices of s and t (in the <p^ 's). These indices could be represented by 0(ln «)-bit 
opaque vertex ID's and could be carried in the message headers. Routing schemes that support opaque 
vertex addressing are called name-independent. Name independence allows for vertex name persistence 
across time (i.e. changing graph topology) and across multiple co-existing routing schemes. 

Computation We use an idealized computational model to facilitate this exposition. The vertices of G 
are viewed as processors, synchronized by a global step counter. During a time step, pairs of processors 
can exchange messages of arbitrary (finite) size as long as they are connected by an edge. We describe 
an algorithm for computing approximations !p^ to all <^" ] in 0(\nn/A) steps, where A is the Fiedler 
eigenvalue of G (the smallest non-zero eigenvalue of L). If G is an expander, then A = 0(1). At every 
step the algorithm sends messages consisting of 0{n) real numbers across every edge and performs 
0(de.g(v) • n) arithmetic operations on each processor v. Using standard techniques, this algorithm can 
be converted into a relatively easy-to-implement asynchronous one. (We omit this detail from here.) It 
is assumed that no graph changes occur during the computation of vertex tables. 

A vector f € R y is distributed if <f„ is stored at v, for all v. A matrix M € K VxV is local (with 
respect to G) if M u , v + implies u ~ v or u = v. It is straightforward that if f is distributed and M is 
local, then M£ can be computed in a single step, resulting in a new distributed vector. Extending this 
technique shows that for any polynomial q(-), the vector q(M)g can be computed in deg(g) steps. 

The Power Method gives us a matrix polynomial q(-) of degree 0(\n n/A) such that q(L) is a "good" 
approximation of iJ. We compute the distributed vectors := q(L)x w , for all w, in parallel. As a 
result, each vertex u obtains <p^ = (Q , . . . ,£u), which approximates tp^ according to Theorem 1.3 
and the symmetry of L. In one last step, every processor u sends </> M to its neighbors. The approximation 
error n~ 5 is chosen to suffice (in accordance with Corollary 5.2) as discussed next. 

THEOREM 1.3 — Let A be the Fiedler (smallest non-zero) eigenvalue of G's Laplacian L, and 
let G be of bounded degree d max . Then \\£ [v] - <p [v] \\ 2 < n' 5 , where £ [v] = (2J max )" 1 £* =0 M">Xv and 
M = I - L/2J max , as long as k ^ Q.(A~ l ■ Inn). 

Robustness and latency In order to get a handle on the analysis of routing in an ever-changing net- 
work we use a simplifying assumption: the graph does not change during the computation phase while 
it can change afterwards, during the routing phase. This assumption is informally justified because the 
computation phase in expander graphs (which we consider to be the typical case) is relatively fast, it 
takes 0(lnn) steps. The routing phase, on the other hand, should be as "long" as possible before we 
have to recompute the scheme. Roughly, a routing scheme can be used until the graph changes so much 
from its shape when the scheme was computed that both the probability of reaching destinations and the 
congestion properties of the scheme deteriorate with respect to the new shape of the graph. We quantify 
the robustness of electric routing against edge removals in the following two theorems: 
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THEOREM 1.4 — Let G be an unweighted graph with Fiedler eigenvalue A = 0(1) and max- 
imum degree J max , and let / [i f] denote the unit electric How between s and t. For any < p < 1, 
let Q p = {e e E : |/ e | > p) be the set of edges carrying more than p How. Then, \Q p \ < 
rmn{2/(Ap 2 ),2dm^\\L%^i/p}. 

Note that part one of this theorem, i.e. \Q p \ < 2/(Ap 2 ), distinguishes electric routing from simple 
schemes like shortest-path routing. The next theorem studies how edge removals affect demands when 
"the entire graph is in use:" 

THEOREM 1.5 — Let graph G be unweighted of bounded degree d max and vertex expansion a. 
Let f be a routing of the uniform multi-commodity demand set over V (single unit of demand between 
every pair of vertices ), produced by an n-competitive oblivious routing scheme. Then, for any < x < 1 , 
removing a x-fraction of edges from G removes at most x ■ (n ■ d max • In n ■ a~ 1 )-fraction of flow from f. 

The expected number of edges traversed between source and sink reflects the latency of a routing. 
We establish (Proven in Appendix D): 

THEOREM 1.6 — The latency of every electric walk on an undirected graph of bounded degree 
^max and vertex expansion a is at most 0(min{/n 1 ' /2 , d max a~ 2 In n}). 



Analysis The main hurdle is Theorem 1.1, which we attack in two steps. First, we show that any 
linear routing scheme R (i.e. scheme for which the operator R : R y — > R £ is linear) has a distinct 
worst-case demand set, known as uniform demands, consisting of a unit demand between the endpoints 
of every edge of G. Combinging this with the formulaic expression for electric flow (1.3) gives us 
an operator-based geometric bound for n^, which in the case of a bounded degree graph is simply 
?7£ < OdlL^li^i) where the operator norm || • ||i^i is defined by HAHi-^i := sup v;i0 ||Ax||i/||jc||i. This is 
shown in Theorem 3.1. Second, we give a rounding-type argument that establishes the desired bound 
on ||L^||i^i. This argument relies on a novel technique we dub concurrent flow cutting and is our key 
technical contribution. This is done in Theorem 4.1. This concludes the analysis of the congestion 
properties of electric flow. 

The computational procedure for the vertex potentials <p^'s (above) only affords us approximate 
versions ip^ with £2 error guarantees. We need to ensure that, when using these in place of the exact 
ones, all properties of the exact electric flow are preserved. For this purpose, it is convenient to view 
the electric flow as a distribution over paths (i.e. the electric walk, defined above) and measure the 
total variation distance between the walks induced by exact and approximate vertex potentials. This is 
achieved in Theorem 5.1 and Corollary 5.2. It is then easy to verify that any two multi-commodity flows, 
whose respective individual flows have sufficiently small variation distance, have essentially identical 
congestion and robustness properties. 
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Related work Two bodies of prior literature concern themselves with oblivious routing. One focuses 
on approximating the shortest-path metric [16, 15, 1, 2], the other focuses on approximating the min- 
imal congestion universally across all possible demand sets [12, 8]. The algorithms in these works 
are essentially best possible in terms of competitive characteristics, however they are not distributed 
and do not address (competitive) performance in the presence of churn. It is not obvious how to pro- 
vide efficient distributed variants for these routing schemes that are additionally resistant to churn. The 
primary reason for this are the algorithmic primitives used. Common techniques are landmark (a.k.a. 
beacon) selection [16, 15], hierarchical tree-decomposition or tree-packings [12]. These approaches 
place disproportiantely larger importance on "root" nodes, which makes the resulting schemes vulner- 
able to individual failures. Furthermore, these algorithms require more than (quasi-)linear time (in the 
centralized model), which translates to prohibitively slow distributed times. 

We are aware of one other work in the theoretical literature by Goyal, et al. [6] that relates to 
efficient and churn-tolerant distributed routing. Motivated by the proliferation of routing schemes for 
trees, they show that expanders are well-approximated by the union of 0(1) spanning trees. However, 
they do not provide a routing scheme, since routing over unions of trees is not understood. 

Concurrently with this paper, Lawler, et al. [10] study just the congestion of electric flow in isolation 
from other considerations like computation, representation or tolerance to churn. Their main result is 
a variant of our graph expansion-based bound on ||L'||i_»i, given by Theorem 4.1. Our approaches, 
however, are different. We use a geometric approach, compared to a less direct probabilistic one. Our 
proof exposes structural information about the electric flow, which makes the fault-tolerance of electric 
routing against edge removal an easy consequence. This is not the case for the proofs found in [10]. 

Organization Section 2 covers definitions and preliminaries. Section 3 relates the congestion of elec- 
tric flow to ||L'||i^i. Section 4 obtains a bound on ||L^||i_^i by introducing the concurrent flow cutting 
method. Section 5 relates the electric walk to the electric flow and proves (i) stability against perturba- 
tion of vertex potentials, (ii) latency bounds and (iii) robustness theorems. Section 6 contains remarks 
and open problems. 

2 Preliminaries 

The Spectral Theory of graphs is comprehensively covered in [3]. The object of interest is a graph 
G = (V, E) (with V = [n] and \E\ = m) undirected, positively edge-weighted by w luu > 0, and not 
necessarily simple. Whenever we use unweighted graphs, we have w uv — 1 if w ~ v and w u v - 
otherwise. The Laplacian is positive semi-definite, L > 0, and thus can be diagonalized as L - UAU*, 
where U € R" x " is unitary and A - diag(/li, . . . , A n ). By convention, we write < A\ < A2 < • • • < A n . 
For every G, A n < 2D, where D is the maximum degree. When G is connected, Ker(L) = 1, and so 
LlJ - UL - 7r ± i where nw denotes projection onto W and iJ denotes the pseudo-inverse of L. On 
occasion we use A^m '■= A2 and A max := A n . The vertex expansion of an unweighted G is defined as 
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min{|5|,|sC||- 



3 The geometry of congestion 

Recall that given a multi-commodity demand, electric routing assigns to each demand the corresponding 
electric flow in G, which we express (1.3) in operator form £(© T c? r ) := WBlJ(® T d T ). Electric routing 
is oblivious, since £(© r <i T ) = © r £(J r ) ensures that individual demands are routed independently from 
each other. The first key step in our analysis, Theorem 3.1, entails bounding by the || • matrix 
norm of a certain natural graph operator on G. This step hinges on the observation that all linear routing 
schemes have an easy-to-express worst-case demand set: 

THEOREM 3.1 — For every undirected, weighted graph G, let fl = W II2 BI3b* W 1/2 , then 



Proof of Theorem 3.1. It is sufficient to consider demand sets that can be routed in G with unit conges- 
tion, since both electric and optimal routing scale linearly with scaling the entire demand set. Let ® T d T 
be any demand set, which can be (optimally) routed in G with unit congestion via the multi-commodity 
flow © T / T . Thus, d T - Y>e fr,eB e , for all t. 

The proof involves two steps: 



Step (i) shows that congestion incurred when routing ® r d T is no more than that incurred when 
routing G's edges, viewed as demands, through G: 



^ £ <liw 1/2 nw- 1/2 i| 1 ^i. 



(3.1) 



||£«Mr)|| G < ||£(e e uvB e )|| G = \\w 



1/2^-1/2 



1->1 





(i) 



e 



e 





J 



J 





J 





T 




T 



\E{® e w e B e )\ G 
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||£(© e «; eJ B e )|| G ( '= 1} \\E(® e w e B e )*W 

(1 = 3) \\wbiSb*ww- 



-i 



-l 





= W 



l/2 uw- l/2 h 



(ii) 



Remark 3.2 — Note that the proof of step (i) uses only the linearity of £ and so it holds for any 
linear routing scheme/?, i.e. one has ||/?(© T <i T )||G < \\R(® e w e B e )\\Q. 

Using Theorem 3.1, the unconditional upper bound in Theorem 1.2 is simply a consequence of 
basic norm inequalities. (See Appendix A for a proof.) Theorem 1.1 provides a much stronger bound 
on ?/£ when the underlying graph has high vertex expansion. The lower bound in Theorem 1.1 is due to 
Hajiaghayi, et al. [7]. They show that every oblivious routing scheme is bound to incur congestion of at 
least Q.(\nn/ In Inn) on a certain family of expander graphs. The upper bound in Theorem 1.1 follows 

from Theorem 3.1, Theorem 4.1 and using that | |TT| 1 1 > i = 0(||L'||i^i) for unweighted bounded-degree 

graphs. Thus in the next section we derive a bound on ||L' in terms of vertex expansion. 



4 L\ operator inequalities 

The main results here are an upper and lower bound on ||Z^||i_>i, which match for bounded-degree 
expander graphs. In this section, we present vertex expansion versions of these bounds that assume 
bounded-degree. 

THEOREM 4.1 — Let graph G - (V,E) be unweigthed, of bounded degree d max , and vertex 
expansion 



The proof of this theorem (given in the next Section) boils down to a structural decomposition of 
unit (s, ?)-electric flows in a graph (not necessarily an expander). We believe that this decomposition is of 
independent interest. In the case of bounded-degree expanders, one can informally say that the electric 
walk corresponding to the electric flow between s and t takes every path with probability exponentially 
inversely proportional to its length. We complement Theorem 4.1 with alower bound on ||L^||i^i proven 
in Appendix B: 

THEOREM 4.2 — Let graph G = (V,E) be unweighted, of bounded degree d max , with metric 
diameter D. Then, ||L^||i_»i ^ 2D/d max and, in particular, ||Z^^||i_»i > (21nw) • (<f max lnc? max ) _1 for all 
bounded-degree, unweighted graphs with vertex expansion a = 0(1). 



a = mm 




then 




(4.1) 
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4. 1 Proof of upper bound on 1 1 U | 1 1 _» i for expanders 

Proof of Theorem 4.1. Reformulation: We start by transforming the problem in a more manageable 
form, 

t \0u\\\ + W n - 1 + 

II^Hi^i := sup = maxUL^Ih < maxHL^,. - Xt )\\u (4.2) 

y+o \\y\\\ w n S*t 

where the latter inequality comes from 

HL^IIi = \\l3ns_txsh = \\n l YlJ(Xs-Xt)h < — max \\l)(x s -Xt)\\i- 

t£ n ' 

Pick any vertices s + t and set if/ = V(Xs ~ Xt)- In light of (4.2) our goal is to bound ||t^||i. We 
think of <A as the vertex potentials corresponding to electric flow with imbalance Xs - Xt- By an easy 
perturbation argument we can assume that no two vertex potentials coincide. 

Index the vertices in [n] by increasing potential as < • • • < i// n . Further, assume that n is even 
and choose a median cq so that xpx < • • • < 1/^/2 < Co < 1/^/2+1 < ••• < iff n . (Jf nis odd, then set cq to 
equal the pottential of the middle vertex.) 

We aim to upper bound ||^||i, given as ||i/HJi = Z„|(A„| = Z K ^„>o^ - E» : ^<o^- Using that 
fiw - 0, we get ||(AHi = 2 £k^>0 </r„ - -2 Ek:^<o *Pu- 

Assume, without loss of generality, that < cq, in which case 

n/2 

Mh = -2 J < 2 £|<fc - col =: 2tf (4.3) 

u:i// u <0 (=1 

In what follows we aim to upper-bound N. 

r 

Flow cutting: Define a collection of cuts (5,-,5 . ) of the form 5, = {v : tA„ < c,}, for integers 
i > 0, where S ,• will be the smaller side of the cut by construction. Let kj be the number of edges cut by 
(Si, S P) and pij be the length of the / h edge across the same cut. The cut points a, for i > 1, are defined 

according to c; - c,_i - A,_i, where A,_i := 2£; ^-—^L. The last cut, (S r+ i,S^, ,), is the first cut in the 

sequence Co, c\, . . . , c r+ i with & r+ i = or, equivalently, S r+ \ = 0. 

Bound on number of cuts: Let = |5/|. The isoperimetric inequality for vertex expansion (4.1) 
applied to (S,,Sp) and the fact that < «/2, by construction, imply 

- > a. (4.4) 

Let /,• be the number of edges crossing (S ,-, 5 . ) that do not extend across c,+i , i.e. edges that are not 
adjacent to S;+i. The choice A; := 2YjjPijlki ensures that /, > fc,/2. These edges are supported on at 
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least li/dmax vertices in 5/, and therefore tii+\ < - /,7<i max . Thus, 

h k i (4 - 4) an i i, a \ 

m+\ < tij - - — < rii - — — < tij - — — - 72,(1 - — — j, (4.5) 

"max -^"max -^"max -^"max 

Combining inequality (4.5) with n$ = n/2, we get 

rii a \i 

m < - 1 - ^- (4.6) 

The stopping condition implies S r + 0, orn r > 1, and together with (4.6) this results in 

r<log 1/e !J, with0= (4.7) 



*max 



Amortization argument: Continuing from (4.3), 

11,1 (*) r ! ' 

N = Y, Wi ~ col < Yj {ni ~ Z A > (4 " 8) 

i=i (=o y=o 

where (*) follows from the fact that for every vertex v € S i - 5, + i we can write \ifr B - cq\ < Z;=o A/- 

Because BlJ(x s - Xt) is a unit flow, we have the crucial (and easy to verify) property that, for all i, 
Yjj Pij = 1- In other words, the total flow of "simulatenous" edges is 1. So, 

v-i Pa 2 (4-4) 2 
A, = 2Y^ = - < — (4.9) 
z y J ki ki ani 

Now we can use this bound on A ; - in (4.8), 

J> - n i+1 ) J] Aj ( = } J] RjAj < § £ 1 = + D. 

r'=0 ;=0 t=l i=0 

where to derive (*) we use n r+ \ - 0. Combining the above inequality with (4.7) concludes the proof. I 



5 Electric walk 

To every unit flow / e R £ , not necessarily electric, we associate a random walk W = Wo, W\,... called 
theyZow walk, defined as follows. Let cr := B*f and so cr v - 0. The walk starts at Wo, with 

tdiw i 2-max{0,cr u } £w /«-»» ~ Z«, fw^v 
1P{W = v} =- 



2jio \°~w\ ( X id ./»—>!« Hiwfw^u) 
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If the walk is currently at W t , the next vertex is chosen according to IP{VKr + i = v \ W t = u) = 



fu—>, 



V 



w Ju 



■w 



where 




\f(u,v)l ( M > v) e E and / flows from u to v 
0, otherwise. 



(5.1) 



When the underlying flow / is an electric flow, i.e. when / = 8.(B*f), the flow walk deserves the special- 
ized name electric walk. We study two aspects of electric walks here: (i) stability against perturbations 
of the vertex potentials, and (ii) robustness against edge removal. 

5.1 Stability 

The set of vertex potential vectors <p^ = L*Xv, for all v e V, encodes all electric flows, as argued 
in (1.4). In an algorithmic setting, only approximations </3 M of these vectors are available. We ensure that 
when these approximations are sufficiently good in an £2 sense, the path probabilities (and congestion 
properties) of electric walks are virtually unchanged. The next theorem is proven in Appendix C: 

THEOREM 5.1 — Letip lv] be an approximation of<p [v] , for all v € G, in the sense that 



where A > 4 is a constant. Then for every electric walk, defined by vertex potentials <p = a v (p^ v ', the 
corresponding "approximate" walk, defined by vertex potentials ip = Yiv a vy> induces a distribution 
over paths y with 



where y ranges over all paths in G, and TP V { W - y) denotes the probability of y under <p (respectively 
forJP^W = y}). 

As shown in Theorem F. 1, the Power Method affords us any sufficiently large exponent A, say 
A = 5, without sacrificing efficiency in terms of distributed computation time. In this case, the following 
corollary asserts that routing with approximate potentials preserves both the congestion properties of the 
exact electrical flow as well as the probability of reaching the sink. 

COROLLARY 5.2 — Under the assumptions of Theorem 5.1 and A = 5, the electric walk defined 
by vertex potentials <p = ip^ - <p^ reaches t with probability 1 - o„(l). Furthermore, for every edge (u, v) 
with non-negligible load, i.e. \<p u - £p v \ = u>(n~ 2 ), we have \£p u - <p v \ — \(p u - <p v \, where <p = <p^ - <p^\ 

5.2 Robustness 

We prove Theorem 1 .4 here, since its proof interestingly relies on the flow cutting techniques developed 
in this paper. Theorem 1.5 is proved in Appendix E. 



\\<p [vi - <p M \\2 < v, for all veV, with v = n 



-A 



(5.2) 




(5.3) 
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Proof of Theorem 1.4. For the first part, let {{u u v{), (u k , v k )} - Q p and let p t - - \(Xm ~ 

X v )*L^{x s - Xt)V Consider the embedding C, : V — > R, defined by = xl^(6 s - 8 t ). Assume for 
convenience that < £(d,-) for all /. Let £ m i n = min„ and £ max = max„ £(i>). 

Choose c uniformly at random from [£ m i n , £max] and let X; = pi ■ I{£(w,-) < c < ^(f,)}, where I{-} is 
the indicator function. Observe that the random variable X = Z,X ; - equals the total electric flow of all 
edges in Q p cut by c. Since these edges are concurrent (in the electric flow) by construction, we have 
X < 1 . On the other hand, 

ex = V Pi • pfc <«,) < c < avd} > Y Pi j p ' >Y^ >k K- 

j ■ bmax hmm ■ * 

Combining this with EX < 1 produces \Q p \ < 2/(Ap 2 ). 

For the second part, kp < ^Q p \f\* A \ < Z ee£ l/i^l - -*t)lll < 2</ max • ||Z.t||i_i. This 

gives \Q p \ < 2<i max • HL 1 "!!^!//?. I 

6 Conclusions 

Our main result in Theorem 1 . 1 attests to the good congestion properties on graphs of bounded degree 
and high vertex expansion, i.e. a = 0(1). A variation on the proof of this theorem establishes a 
similar bound on t]£, however, independent of the degree bound and as a function of the edge expansion 

P = mins = y minlS).'voi(5C)) ' where vol(5) : = ^veS w u,u and vol(E(S,S^)) := Z^^^sC «W- 

The bounded degree assumption is also implicit in our computational procedure in that all vertices 
must know an upper bound on d max in order to apply M in Theorem 1.3. Using a generous bound, 
anything o»(l), on d max is bad because it slows down the mixing of the power polynomial. To avoid this 
complication, one must use a symmetrization trick outlined in Appendix G. 

We conclude with a couple of open questions. A central concern, widely-studied in social-networks, 
are Sybil Attacks [4]. These can be modeled as graph-theoretic noise, as defined in [9]. It is interesting 
to understand how such noise affects electric routing. We suspect that any 0(ln ^-competitive oblivious 
routing scheme, which outputs its routes in the "next hop" model, must maintain £2(7i)-size routing tables 
at every vertex. In the next hop model, every vertex v must be able to answer the question "What is the 
flow of the (s, f)-route in the neighborhood of vT in time C?(polylog(?i)), using its own routing table 
alone and for every source-sink pair (s, t). 
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A Proof of Theorem 1.2 

Proof of Theorem 1.2. The upper bound follows from: 

(3-D ( A - 2 ) 1/9 (A3) ,/9 

77 < IIITNi^i < m 1/2 • ||H||2— 2 = m 112 (Al) 

The second step is justified as follows 

||n||i^i = max ||ILy e ||i < m 1/2 • max ||n^|| 2 < m 1/2 • ||n|| 2 ^ 2 . (A.2) 

e e 

The third step is the assertion 

linn*-* < i, (A.3) 

which follows from the (easy) fact that n is a projection, shown by Spielman, et al. in Lemma A.l. 

The lower bound is achieved by a graph obtained by gluing the endpoints of yjn copies of a path 
of length yfn and a single edge. Routing a flow of value -\/n between these endpoints incurs congestion 
y/n/2. I 

LEMMA A.l (Proven in [14]) — II is a projection; Im(II) = lm(W 1,2 B); The eigenvalues of IT 
are 1 with multiplicity n - 1 and with multiplicity m - n + 1; and U e , e = ||riY e || 2 . 

B Proof of Theorem 4.2 

Proof of Theorem 4.2. Let s t t be a pair of vertices in G at distance D. We consider the flow / - 
BlJ(Xs - Xt)- Set tA = L^ixs - Xt), and note that we can use ||/||i as a lower bound on ||L^||i^i, 

\\f\h = £ \*u ~ <A»I < d^^vl < d max \\L^ s - Xt )\h < ^II^Hi^i. 

(u,v) V ~ 

Now, let be a path decomposition of / and let l(m) and /(tt/) denote the length and value, respec- 
tively, of 7T(. Then, 

\\f\U = 2 Wu - 4>v\ = K*i)f(fii) f{m) = D. I 

(u,v) i i 
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C Proof of Theorem 5.1 

Proof of Theorem 5.1. Notation: Note that in this proof we use the notation of (5. 1). So, for a potential 
vector iff, we have \p u ~, v := \f/ u - iff v if (u,v) is an edge and if/ u ^ if/ v , and ifr u -> v := otherwise. So, 
for example, the potential difference on (u, v) can be written as iff u -+ v + ifr v -+ u - On the other hand, we 
use the single letter edge notation <p e to denote the signed (according to B) potential difference on e, so 
(p e := ip u - ip v if B e = 6 U - S v . Let D be the maximum degree. 

Edge approximation: Fix any unit electric flow, defined by potentials tp :— Y* v a vf^ v \ and write its 
approximation as £p := a v ft s \ All unit flows can be so expressed under the restriction that a v = 
and Yj V \ a v\ - 2- The approximation condition (5.2) combined with Lemma C.l gives us, for every edge 
e = (u, v), 

\<Pe ~ <fie\ = | 2j a v(f^ ] ~ ^e" ] )| 

V 

,r, i i m ~mi 

< ^ |ar u | • 2y apply Lemma C. 1 

V 

= 4v 

We call this the additive edge approximation condition 

(fie - 4v < ip e < (f e + 4v (C.l) 

Now, consider a fixed path y along the electric flow defined by tp, traversing vertices wq, w\, . . . , 
Let P V {W = y) and P^{VK = y} denote the probability of this path under the potentials tp and tp, 
respectively. In most of what follows, we build machinery to relate one to the other. 

Path probabilities: For a general unit flow (not necessarily an (s, 0-flow), defined by vertex poten- 
tials \fr, IP|/,{W = y] equals 

EVIWo = wJ Y\ V*{W i+1 = w M I Wt = w^Y^Woo = w k | w k }, (C.2) 

where next we explain each factor in turn. 

The first, IP^{Wo = wo), is the probability that the walk starts from wo, and is expressed as 

IV { Wo - w \ = max (0, ^ <p wo -> u - 2_j fu-^w)- (C.3) 

u u 

The second and trickiest, P^{W,-+i = Wj+i \ Wi - Wi), is the probability that having reached Wj the 
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walk traverses the edge leading to Wi + \, and £ M iffuii-m ^ E« ( Au^to,> we write 

P^{W/ + i - 1 Wi - a;,-} ^"* + ' ^_ (C 4) 

max ( 2^ tu^, 2j ^wr+u) 

U U 

To grasp the meaning of the denominator, note that the quantity 1 2a ftu—*wi ~ 2j« tywi— m\ is the magnitude 
of the in or out flow (depending on the case) at &>,-. 

The third, IP^IWoo = Wk I Wk\, is the probability that the walk ends (or exits) at wu conditioned on 
having reached wu, and 

IP^{ Woo = w k I w k ) = max (0, ^ ifr u -+w k ~ 2_j ^«*-«)' (C.5) 

Next, we are going to find multiplicative bounds for all three factors by focusing on "dominant" 
paths, and discarding ones with overall negligible probability. 

Dominant paths: It is straightforward to verify (from first principles) that the probability that an 
edge (u, v) occurs in the electric walk equals \<p e \ = + <p v ^> u - We call an edge short if \<p e \ < e, 
where the exact asymptotic of e > is determined later, but for the moment v «: e «: 1 . We restrict our 
attention to dominant paths y that traverse no short edges, and have fl^Wo = wq\ ^ £ and IP^VKoo = 
Wk I w k ) > e. 

Indeed, by a union bound, the probability that the electric walk traverses a non-dominant path is at 
most 2ne + n 2 e. This will be negligible and such paths will be of no interest. In summary, 

JP^W dominant} > 1 - 2ne - n 2 e (C.6) 

We now condition on the event that y is dominant. 

The no short edge condition gives e < \tp e \ < 1, and using (C. 1) we derive the stronger multiplicative 
edge approximation condition 

1 w e 8Dv 

- < — < cr, where cr - 1 + , (C.7) 

cr <fie e 

which holds as long as e ^ 4v, as guaranteed by the asymptotics of e. Also note that the latter condition 
ensures that tp e and ip e have the same sign. An extra factor of 2D is included in cr with foresight. 

For the first factor (C.3), we have 

JP<p{W = Wq] = ^ (fw ^u ~ ^ <Pu-*w (C8) 

U U 

> 2j Vm^u ~ ¥i^w ~ 4Dv use (C. 1) 

u u 

> TP v {Wo = w }(l - — ) use V^Wq - w ] > e 

>-JPJW = w Q } use e< 1/2. 

cr 
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For the second factor (C.4), assume 'Z u i// u ->w i > H« <A»;,— >«• An identical argument holds in the 
other case. Abbreviate 

Pp{u>i+i | := Hy - w M | Wi = wtl 
Path dominance implies Y* u fu-^wt ~> £, and so 

Py{«;,- + 1 1 «;,-}= ^ Wi+l (C.9) 



; (Pwj-+w i+l 



> — use (C.7) and (C.l) 



> cr = use 6 < 1/2 and ^fu^w, > e 



For the third factor (C.5), similarly to the first, we have 

P^{Woo=lDfc|lD t }= (CIO) 
u u 

> 2_j <pu-*w k - 2] ^w*^" - 4£)y use ( c - !) 

U W 

/ 4Dvx 

5* P v {Woo = | w k }[\ —J use P y {Woo = Wk\wk)> e 

> —IPJWoo = w k \w k } use e< 1/2. 

cr 

Dominant path bound: We now obtain a relation between P v { W - y) and P^{W = y) by combing- 
ing the bounds (C.8), (C.9) and (CIO) with (C.2): 

Py{W = y] 1 

^— tttt r > -^rp? a PP J y bounds ' and P ath len § th < " (CM) 

I 8Dv\2n+2 , 

> (l ) use cr -1 ^ 1 - 8Dv/e 

/ 16Dv\2n+2 _ 9r 

> exp y j use 1 - x > e 



Statistical difference: Abbreviate p(y) := P V {W - y} and q(y) := P V {W - y}. Below, y iterates 
through all paths, f iterates through dominant paths and £, iterates through non-dominant paths. We 
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bound the statistical difference (5.3), using (C.l 1) which says q{£) > 6 ■ p{£), 

2 \p{y) - -7(7)1 - 

7 

= £ \p(0 - q(0\ + £ \P(0 ~ q(D\ 

< 2 1(1 - 0M<T) - (q(0 - 9p(0)\ + £ q{& use < 1 



In this final step, we pin-point the asymptotics of e that simultaneously minimize the two terms of (C. 12). 
In the following, we parameterize e = n~ B and use (C.6), 



(1-0) + Q^PiD 



, 16Dvx2n+2 / 16Dvv2n+2 , 

= l-exp^ j + exp^ j (2ne + n e) 

= 1 - exp 0( - Dn B ' A+1 ) + n 2 ' B ■ exp 0( - Dn B ~ A+1 ) 

= 0{Dn B ~ A+l ) + n 2 ~ B ■ exp 0( - Dn B ~ A+l ), use 1 - e~ x < x 

- 0(n B ~ A+2 ) + 0(n 2 ' B ) use D < n 

= <9(ra 2 ~5), set 5 = A/2. 
LEMMA C.l — Ifx,ye( 2 and \\x - y\\ 2 < v, then for all i ± j, 

(x { - Xj) - 2v < y\ - yj < (Xj - xj) + 2v. 



Proof. We have (x ; - - yj) 2 < \\x - y\\ 2 < v 2 , implying |x, - yi\ < v. Similarly for j. Combining the two 
proves the lemma. I 
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D Proof of Theorem 1.6 

Proof of Theorem 1.6. Let x\ s ''^ be the indicator that edge e participates in the electric walk between s 
and t. Then the latency can be expressed as 

max ^ EX [ e SJ] = max ^ \&u - S^)l)(5 s - 6,)\ 

e (u,v) 

= wax\\BL\5 s -5 t )h = H^li^i - Mi-n- 

The latter is bounded by Theorem 4.1 and ||n||i^i < m 1 / 2 , as in (A.l) e.g. I 

Remark D.l — For expanders, this theorem is not trivial. In fact, there exist path realizations of 
the electric walk which can traverse up to 0(n) edges. Theorem 1.6 asserts that this happens with small 
probability. On the other hand, in a bounded-degree expander, even if s and t are adjacent the walk will 
still take a 0(log «)-length path with constant probability. 



E Proof of Theorem 1.5 



Proof of Theorem 1.5. Let / Q pt be a max-flow routing of the uniform demands and let 9 be the fraction of 
the demand set that is routed by / opt . The Multi-commodity Min-cut Max-flow Gap Theorem (Theorem 
2, in [11]) asserts 

nn ^ a ^ ■ g^j^ ^ 1 ■ l £ (^ C )l « 
0(m n) • 9> mm — > - • mm — = — 

scv \S\.\SC\ n scv m i n {\ S \,\sZ\) n 

Thus the total demand flown by / opt is no less than 9^ > Q(an/lnn). Normalize / (by scaling) so it 
routes the same demands as / opt . If k edges are removed, then at most nk flow is removed from /, which 
is at most a fraction rjk ■ 0(lnn/an) of the total flow. Substitute x = k/m and use m < d max n to complete 
the proof. I 



F Proof of Theorem 1.3 

Theorem 1.3 is implied by the following theorem by specializing e = 0(n~ 5 ): 

THEOREM F. 1 — Let G be a graph, whose Laplacian L has smallest eigenvalue A and whose 
maximum degree is D. Then, for every y with \\yW2 = 1 the vector x = l)y can be approximated using 

1=0 
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so that for every e > 0, 



\\x - x\\2 < e, as long asd> • l n 



AeD v I - A' 



Proof of Theorem F.l. We normalize L via N = L/t (and so L 1 = N 1 /t), where r = 2D. Since 
t = 2D ^ /l max (L), the eigenvalues of N are in [0, 1]. In this case, the Moore-Penrose inverse of N is 
given by N~ l = X" ( 7 ~ Set = ^f=o( 7 ~ A i and M = ^ ~ 0ur aim is to minimize so that 



II* - X\\2 



nI + n\ 



Nl 



-y y 

T T 



N, 



2^2 



where ||A||2->2 := su Px^o IIAxlh/IWb denotes the matrix spectral norm. Set k :- r/A m m, so that k 1 is the 
smallest eigenvalue of N, 



l|tff|| 2 -2 = II E (/_ N)i Wz-*2 < Z 11(7 -^112-2 
i=d+l i=d+l 

CO 

< Z (l-^Hl-*- 1 )^ 1 * 

Setting (F. 1) less than re gives 



i=d+l 



(F.l) 



d ^ 



ln k/(t€) 
1iuc/(k- 1)' 



G Symmetrized algorithm 

In this section we discuss how to modify the computational procedure, given in the Section 1 , in order 
to apply it to graphs of unbounded degree. The described algorithm for computing tp^ = l^Xw relies 
on the approximation of l) via the Taylor series yz^ = Z^o^ _ x )'- The sef i es converges only when 
1Mb < 1, which is ensured by setting x = and using that ||L||2->2 < 2<i max . Thus we arrive at 

^"max 

2d max ■ iJ = Z~o( 7 _ iT~ ^his a PP roacri continues to work if we replace d max with any upper bound 
hmax ^ d max , obtaining L) = t^J— T^lM - MY where M - however this is done at the expense of 
slower convergence of the series. Since in a distributed setting all vertices must agree on what M is, a 
worst-case upper bound h max = n must be used, which results in a prohibitively slow convergence even 
for expander graphs. 

Instead, we pursue a differnt route. Let L = D~ l/2 LD~ { t 2 be the normalized Laplacian of G, where 
D € R nx ' 1 is diagonal with D v v - deg(u). One always has ||£||2->2 < 2 (Lemma 1.7 in [3]) while at the 

same time A m { n {L) ^ max (%-, — } (Theorems 2.2 and 2.6 in [3]), where a and B are the vertex- 

and edge-expansion of G, respectively. Set M = XL/3, so that ||M||2-»2 < 1- Recall that the aim of our 
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distributed procedure is to compute cp\^ at u (for all w). We achieve this using the following: 

Yu - Xu 1 ^ Xw - Xu u o u Xm - r — o r — . 

* v de g( M ) J V de g(" ; ) 

The key facts about the series in the left-hand side are that (i) it converges quickly when G is an expander 
and (ii) all vertices can compute M locally, in particular, without requiring any global knowledge like 
e.g. an upper bound on <i max . 
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